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Similarity renormalization group procedure identifies the role of bound states in the low-energy 
rate of change of effective coupling constant in a model Hamiltonian with asymptotic freedom. 



I. INTRODUCTION 



The fact that an asymptotically free interaction blows up at low energies is often associated with the concept of 
confinement. It is therefore of interest to construct a soluble model that exhibits a similar variation of a coupling 
00 ' constant and see what mechanism causes that the coupling constant rises at low energies. 

In a very simple model with asymptotic freedom and a bound state, similarity renormalization group (SRG) 
procedure for Hamiltonians produces a strong effective interaction at low energies. The interaction strength is related 



o 



to the bound-state formation and a degree to which the interactions that are responsible for the binding are included 
in the generator of SRG transformations: the more interaction in the generator the less increase in the coupling 
' constant Similar results are found in the case with limit cycle, of which the apparent asymptotically free SRG 
behavior of effective interactions may be a part. This lecture concerns only some aspects of these results. Readers 
interested in a bit broader picture may consult Ref. 0. A connection with AdS/CFT correspondence 0, Q that 
inspires a new generation of models for calculating masses of bound states of asymptotically free quarks and gluons 
[1, El, is not yet established. 

O II. MODEL 

<D ' 
(— ! 

Consider a model Hamiltonian in the generic form of 

H = Ho + Hl (1) 

> , 

Ho can be just a free (kinetic) energy or include also dynamical effects that are well understood. The origin of division 
of H into Ho and Hi is that all eigenvalues and eigenstates of Ho are known. Thus, Ho provides a basis for studying 
new effects due to the postulated Hj. For simplicity of the model, suppose the eigenvalue problem for Ho, 
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Ho\n) = E n \n), (2) 



is solved with a set of discrete eigenvalues E n > that form a monotonic sequence, E m > E n when m > n, there is 
no degeneracy, and the corresponding eigenvectors are normalized, (m\n) = S mn . The interaction Hamiltonian Hj is 
^ . then defined by its matrix elements Hj mn — (m\Hj\n). An analytically soluble model to be discussed here, is obtained 
by assuming that these matrix elements have a factorized form, i.e., Hj mn = H m H n . H n should have dimension of 
square root of energy and the simplest choice that does not introduce any dimensionful scale is H n ~ y/E n . The 
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■ proportionality is reduced to a dimensionless number, and one can write 

{m\H I \n) = -g^E~^EZ, (3) 

where g determines the strength of the interaction, called coupling constant in analogy with QFT. The negative sign 
results in existence of a bound state for sufficiently large positive g. By definition, the bound state corresponds to 
a negative eigenvalue of H . It will have to be clarified what the words "sufficiently large g" are supposed to mean, 
because the model Hamiltonian produces divergences (infinities) no matter how small the number g is. 
In summary, the model Hamiltonian is defined by its matrix elements as follows: 

H mn = (m\H\n) = E m 5 mn - g\[E^ l \[E^ . (4) 
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III. ULTRAVIOLET DIVERGENCES 



Suppose g is a very small number and one expects that eigenvalues of H should be nearly equal to the eigenvalues 
of H . First-order correction to E m is a fraction g of the energy being corrected, 



A£« = (m\Hi\m) = -gE m . (5) 



The second-order correction is 



r~r &m — &k f~r — &k 

k^m kytm 

Terms with k < m involve Ek < E m and if Ej~ < E m /b m ~ k with a number b > 1, they contribute an amount that is 
not sensitive to the lower bound on Ek, say b M with a large negative integer Af, when M — ► — oo. Terms with fc 3> m 
involve ^> _E m and thus each of them contributes —1. If Ek = b k ~ m E m , and b 3> 1, the second-order correction 
to is proportional to the number of basis states with energies (eigenvalues of H ) greater than E m . In order to 
obtain a finite correction, one has to limit the number of the basis states and try to understand what happens when 
the limit is relaxed. In analogy with QFT, imposing such a limit is called regularization. If one assumes a cutoff of 
the form k < N, where N is a large positive integer, the result of Eq. © for 6 ^ 1 is 

AE$ = -g 2 (N-m)E m . (7) 

The ultraviolet cutoff on energies, A = b N ~ m °E mo means k < N — mo + \nA/E mo , and one says that the 
second-order correction is ultraviolet divergent logarithmically, 

AE - ^ -^^^e\ e - ^ 

The divergence results from all different energy scales contributing equally to the correction. The occurrence of 
divergences is not limited to perturbation theory. 

With the factorized interaction, solutions to the eigenvalue problem 

N 



^ H mn i/j n = Eip m , (9) 



n=M 

have wave functions 

N 

n 

E m — E 



^ m = V, E -K 9 51 V^An , (10) 



in which the eigenvalue satisfies the condition 



N 



n=M 

The sum here resembles closely the one in Eq. ^ and develops the same type of divergence. If a bound state exists, 
with a negative eigenvalue E = —Eb, one can replace the sum over n by an integral, dn = dE n j {E n In 6), to estimate 
what happens, and obtains 

A-b M e-ir 

Eb - ■ (12) 

e s — 1 

In the limit of A — * oo for fixed g, the binding energy diverges. One can also observe that the square of the matrix 
Hi with large N and M is equal to Hi, which means that it diverges in the ultraviolet limit of A — > oo for 
fixed g. This means that all powers of the entire H are ultraviolet divergent. In particular, the evolution operator 
U(t, 0) = e~ lHt does not exist in this limit. 
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IV. ASYMPTOTIC FREEDOM 



The trouble with divergences can be summarized using the model as follows. One measures transition rates between, 
say, two states, say \rrii) and \rri2), and discovers that these rates can be described in first-order perturbation theory 
by an interaction Hamiltonian with matrix elements 



(m 2 \Hi\m 2 ) 
(mi\H I \m 2 ) 



(m 2 \H I \mi) 
(mAHAmi) 



-9 



(2) 



(13) 



The coupling constant corresponds to physics of the 2 states. One is then compelled to postulate that the whole 
matrix of Hi has the form given in Eq. ([3]) . This step is analogous to the proposal of non-Abelian gauge theory [t| . 
Even if this leads to divergences, one does not want to abandon the proposed interaction since it does produce 
a structure in perturbation theory that fits the case of states \m\) and \rri2) and it has an appealing symmetry. 
Therefore, one looks for a general way out of the problem with divergences that are produced by naive extrapolation 
of knowledge from a small set of matrix elements to a large set. The large set is desired when one seeks a theory of 
presumably large range of applicability and a lot of predictive power. 

The way to proceed is to learn what happens when one begins with some large N and tries to mathematically reduce 
N to a small value near m 2 . As a principle, such procedure was proposed and developed in seminal Refs. [8j, |9fl. I t 
is sometimes called "integrating out high-energy degrees of freedom," which applies also in statistical mechanics 10]. 
In the model, one starts with the eigenvalue problem 



= E[ 



(14) 



k=M 



and applies Gaussian elimination, beginning with ij)ff. The remaining set of equations for ipk with k < N — 1 
corresponds to (the required algebra is merely solving one linear equation) 



H 



(N-l) 



[( N-1) 

the bracket). A simplification occurs for cutoffs much larger than \E\, for which E/Em can be neglected, 



HV 1 ' has the same structure of matrix elements as Hi but contains a new "coupling constant" (the expression in 
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(N-l) 



9 



1 - E/E. 



N 
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1-g- E/E N 1 - g 
Therefore, one obtains a recursion that does not depend on the eigenvalue E, 
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(K-l) 



9' 



l-g(K) 



for as long as Ak = b H> \E\. This is the RG recursion in the model. It is solved by 
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{K) 



(16) 



(17) 



(18) 



l-g{N-K) 

Suppose the eigenvalue problem for H^ K °^ with b K ° — Ao is small enough to solve for its spectrum precisely using 
computers and establish that the coupling constant g( K °> should have some value go in order to reproduce some 
measured eigenvalue Eo with |2?o| *C Ao (some transition amplitude could be used instead). For E n — 6", Eq. (fl8|) 
says that 



5o = 



9A 



1-^lnA/Ao 



A = 



9o 



1 + ^lnA/Ao 



(19) 



This means that the model is asymptotically free [ll|, [l2| : the larger the cutoff A in the initial H the smaller the 
coupling constant <?a in it, and the smaller value is required to describe high-energy processes in perturbation theory. 
Calculation of g\ means also evaluation of the counterterm required in the initial Hamiltonian. 

The increase of g^ when K decreases causes a major difficulty because the ratio E/ Ek is compared with 1 — g( K > 
in the RG recursion and the unknown eigenvalue E cannot be ignored when g( K > — > 1 no matter how much smaller 
\E\ is than Ek- This is a basic obstacle to description of bound states in theories with asymptotic freedom: the 
eigenvalues cannot be found in perturbation theory and one cannot easily reduce the cutoff to sufficiently low values 
for carrying out all kinds of interesting non-perturbative calculations using available computers. In order to control 
what happens in the range where g^ K > is order 1, a different RG procedure seems appropriate. 
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SIMILARITY RG PROCEDURE AND gx AT LOW ENERGIES 



In the SRG procedure [13| , one proceeds according to similar principles as in the standard approach described in the 
previous section. One also finds counterterms and evaluates effective Hamiltonians. The new idea is that one does not 
"integrate out" any degrees of freedom. Instead, one changes the basis states by rotating them in the Hilbert space. 
The rotation is designed in such a way that it guarantees the resulting Hamiltonian matrix, H\, to have vanishing 
matrix elements between basis states if they differ in energy by more than A. The design is such that one does not 
encounter small energy denominators even when the coupling constant increases to 1 or larger values. Of course, the 
SRG procedure produces the same CT as the one identified in the previous section and the initial condition at A = oo 
for the SRG evolution of H\ with A is the same 



H Ar 



— E m 5 mn — g\ 1 



The SRG evolution can be obtained from (prime denotes differentiation with respect to A) 

H'\ = [T\,H\], 



(20) 



(21) 



where T A = [G\,H\], G\ = fH + (1 — f )Dx, and D\ denotes the diagonal part of H\. For / = 0, one has G\ = D\ 7 
in which the diagon al p art of interactions is fully included and Eq. (|2"Tj) is the one introduced by Wegner in condensed 
matter physics [14, KLa]. For / = 1, one has Gx — Hq, in which no interaction effects are included, and Eq. (|2~lj) is 
then the one used in nuclear physics [HI, [13, E3] . For intermediate values of / 6 [0, 1], G\ includes interactions to an 
intermediate degree, correspondingly, and one can inspect what happens in various cases. 
Numerical calculations produce results that can be summarized by writing 



-[E m —E n )/X 



(22) 



The effective coupling constant g\ is then defined using the interaction matrix elements between the states of lowest 
energies, such as |irii) and \7ri2) in Eq. (|13p . in analogy with the Thomson limit in QED, say mi = M , m 2 = M + 1, 
and 



.9a — 1 — H\mm+i/ v EmEm+i ■ 
A generic example is shown in Fig. [TJ where 6 = 4, N = 16, M = —25, g\ = 



9A 



(23) 

4/100, and A = 4 16 - 4 • 10 9 . 




FIG. 1: gx grows toward small A differently for / = 0, 0.5, 0.75, 0.9, and 1 (the larger /, the higher curve). The huge 
increase of g\ below In A/ Info ~ —8 for / = 1, occurs because A decreases there below the scale of binding energy Eb- In fact, 
gx — > \EB\/b M ~ 10 10 . For / = 0, the huge increase of gx is absent. 

The bound-state energy is Eg ~ —8 • 10~ 6 . Fig. [1] demonstrates that the increase of g\ at small energies is caused by 
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removing interactions responsible for existence of a bound state from G\ in the generator of the SRG transformations. 
When the generator fully accounts for the interactions, the magnitude of the coupling constant g\ never significantly 
exceeds 1. 

An apparently very small alteration of Hi mn in the model, by a term —i /isgn(m — n)\/ E m E n with a very small 
coupling constant h, leads to a new way of thinking about asymptotic freedom as a part of a limit cycle. RG limit 
cycles were introduced in the context of strong interactions in Ref. and recently suggested relevant to the infrared 
behavior of QCD [20]. If h is an irrational number, the model typically exhibits a chaotic RG behavior. When 
h = tan - with p a large integer, a limit cycle occurs, with a period \P ~ e*l h . This means that g\ x has the same 

value as g\ 2 if Ai = (b p ) k \2 with integer k. The cycle is associated with existence of bound states whose binding 
energies form a geometric series with quotient l/b p . If a RG cycle were indeed present in a realistic extension of the 
standard model with some tiny coupling constant like h, say k, a new generation of particle substructure would be 
predicted with binding energies order g^^Kqcd- 

VI. CONCLUSION 

The simple model shows that the SRG procedure may be a suitable tool to handle the increase of the coupling 
constant g\ in QCD when A — * Aqcd- If the generator of SRG transformations does not include interactions in G\, 
the effective coupling constant in the model quickly increases to very large values as soon as the SRG scale parameter 
A is lowered down to the momentum scale that characterizes formation of a bound state. This scale is much larger 
than the scale associated with confinement, which concerns distances far beyond the size of a single hadron. If the 
generator of SRG transformations includes interactions in G\, the SRG parameter A can be brought down to the 
momentum scale that characterizes the bound state and g\ does not increase to large values. The SRG procedure 
also helps us recognize a connection between asymptotic freedom and limit cycle, and the model example shows that 
in order to handle the case of limit cycle the generator of SRG transformations must include interactions in G\ . 
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